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Abstract 

In this paper we show how to split the root-locus plot for an irreducible rational transfer function 
into several individual algebraic plane curves, like lines, circles, conics, etc. To achieve this goal we use 
results of a previous paper of the author to represent the Root Locus as an algebraic variety generated 
by an ideal over a polynomial ring, and whose primary decomposion allow us to isolate the planes curves 
that composes the Root Locus. As a by-product, using the concept of duality in projective algebraic 
geometry, we show how to obtain the dual curve of each plane curve that composes the Root Locus and 
unite them to obtain what we denominate the “Algebraic Dual Root Locus”. 

Index terms — Root-Locus, Projective Root-Locus, Ideal of Polynomials, Primary Decomposition, 
Real Projective Plane, Dual Algebraic Curve, Grobner Basis. 


1 Introduction 

Root-Locus (RL) is a parametric plot of the roots of the polynomial p(s) = d(s) + kn(s) over the complex 
plane, equivalently over the affine plane R 2 , as the parameter k spans R; d and n are fixed coprime polynomi¬ 
als, and d is monic with degree, in general, greater than the degree of n. The polynomial p can represent the 
denominator of a (proportional) control feedback loop of a linear time invariant plant with transfer function 
G(s) = n(s)/d(s) (see Figure [lj, and this makes the RL a classical approach to study stability and perfor¬ 
mance of closed loop feedback systems. The rules for sketching the plot are discussed in most textbooks on 
feedback control theory of linear systems (see &)■ In a previous paper ([!]) the author showed that the RL 
plot can be extended to the real projective plane (RP 2 ) and be interpreted as a projective algebraic variety, 
that we denominated projective root-locus (PjRL). In this approach, the RL points, including the ones at 
infinity, can be calculated as the solution of a set of homogeneous polynomial equations, as well as, we can 
obtain complementary plots of RL over the different affine planes that make up the projective plane. 

In this paper we use the approach presented in [I] to show how to decompose the RL into several planes 
curves that can be plotted independently to form the final RL plot. In fact, at least in some simple cases, we 
can easily visualize the RL as a union of several plane curves: for example, the plot presented in Figure [2] 
that represents the RL for G(s) = (s + l)/s 2 , is composed by the (parametrized) circle (x + l) 2 + y 2 = 1 
and by the (parametrized) line y = 0. In order to deal with this question in a systematic way, however, we 
need concepts from algebraic geometry, considering the RL as an algebraic variety, and the goal is to find its 
decomposition into irreducible components (see [BJ Chap. 4]). We cannot, in general, obtain the irreducibles 
components of an algebraic variety “by hand”, but considering the curve as the set of zeros of an ideal in a 
polynomial ring, the question becomes strongly related to the computation of primary decomposition of ideals 
(see JO Chaps. 4,7]), a fundamental topic in abstract algebra, and for which computing algorithms there 



Figure 1: Control Feedback Loop with a Proportional Controller 
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Figure 2: Root-Locus for G(s) = (s + l)/s 2 


exists since a long time ([3]). In particular, Macaulay2 package (Hj) incorporates a command to compute 
the primary decompositon of an ideal. 

We also present in this paper, mainly as a matter of mathematical curiosity, a new root-locus plot that 
we denominate “Algebraic Dual Root-Locus” or (ADRL), associated to the conventional RL plot, that is 
obtained by computing the dual curve (in projective geometry sense) for each individual plane curve that 
makes up the projective root-locus. We leave the analysis of the properties of ADRL for a possible future 
work. 

Bellow we present some concepts used in the paper: 

R, C and R[xi,X 2 , ... ,x n ]: Represents the field of real numbers, the held of complex numbers and the ring 
of polynomials with coefficient’s in R and with indeterminates (xi,X 2 , ■ ■ ■ , x n ), respectively. 

Homogeneous polynomial A polynomial (in several variables) is homogeneous when all of its nonzero 
terms (monomials) have the same total degree (sum of the degree of each variable). We always can 
turn a non-homogeneous polynomial ( q ) into a homogeneous one (q h ) by adding a new variable (x n +i), 
with the following procedure: q h (x 1 ,.. .,x n ,x n+1 ) = x^ +1 q(xi/x n+1 ,x 2 /x n+1 ,.. .,x n /x n+1 ), where 
d is the total degree of q; this process is denominated “homogenization” of q. We can always “de- 
homogenize” q h by setting x n+ \ = 1 and recover back q. 

Ideal of Polynomials: A set of polynomials I C R[xi, X 2 ,..., x n ] is an ideal when it satisfies the following 
properties m, H) : w 0€ /; (ii) p,q £ I implies p + q £ I; and (in) p £ I and q £ R[xi, X 2 ,..., x n ] 
implies pq £ I. One important fact about ideals of the ring R[xi,X 2 , ..., x n ] is that they are finitely 
generated , that is, for every ideal I always there exists a finite subset of polynomials in /, denoted by 
{pi,P 2 , ■ ■ ■,Pt }, such that 

t 

I = ^2hiPi, hi £ R[xi,x 2 , ... ,x n \. 

i=1 

The set {pi,p 2 , ■ ■ ■ ,Pt} is denominated a generating set for /; in this case we write I = (pi,p 2 , ■ ■ ■ ,Pt)- 
A Grobner Basis for the ideal I is a particular kind of generating set that allows many important 
properties of the ideal to be deduced easily. Given a generating set {pi,P 2 , ■ ■ ■ ,Pt} for /, we can 
obtain a Grobner basis {gi, g 2 , ■ ■ ■,g s } for / algorithmically (see 0 Ch. 2]). The most basic ideals 
are (0), the zero ideal, and (1), the ring R[xi, X 2 ,..., x n \ itself. In fact, if 1 £ J, we immediatelly 
conclude I = R[xi,X 2 , ■ ■ •, x„]. We also have that the intersection of any family of ideals results in 
an ideal. Another important result related to ideals in Noetlrerian rings (like R[xi,X 2 , ...,x„]) is the 
“Lask-Noether Theorem” which states that every ideal in a Noetherian ring can be written as an finite 
intersection of primary ideals. 

Variety generated by an ideal: An (real) algebraic variety is a subset of R” whose elements are the (real) 
solutions a system of polynomial equations in n variables (in R[(xi, x’ 2 ,..., x n ]). We can see this set of 
polynomials as a generating set of an ideal J, and so we say that the variety is generated by the ideal 
I, and represented by V(J). There exists several important relationships between ideals and verieties, 
in particular: V((0)) = R n , V((l)) = 0 and if / and J are ideals we have that ¥(/ D J) = V(J) U V(J). 

For more details about the concepts above see (®, 0). 
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2 Decomposing projective root-locus into irreducible components 

In a previous paper (HI) the author showed how to extended the RL plot from the affine plane R 2 to the 
projective plane PM 2 by considering the parametric plot roots of the “modified” polynomial 

p(s) = k d d(s) + k n n(s) (1) 

over PR 2 as k n /kd spans the projective line PR 1 . Considering the ideal I = {u,v), generated by the 
polynomials u = Re{p(a; + iy)} and v = lm\p(x + iy)}, the projective root-locus (PjRL) is obtained from 
the Grobner basis { 31 , 52 , ■ • ■ ,3s} for the ideal /, defined in the ring R [x,y,k d ,k n ], with respect a graded 
monomial order. If we define the homogenization of the ideal I = (u,v) as the ideal I h = ( 31 , 32 , • ■ • ,3s), 
where 3 } is the (homogeneous) polynomial obtained by the homogenization of 3 ^, we can obtain the projective 
root-locus from the (projective) variety generated by the ideal I h , denoted by Y(I h ), where I h is defined in 
the ring R[.t, y, z, kd, k n ] (see [T] for details). 

To decompose the PjRL into irreducible components we need to obtain a primary decomposition of I h . 
in order to write it as a finite intersection of ideals, that is 

i h = j x r\J 2 n ■••n J m , (2) 

where each is a primary ideal (see [TJ Thm. 7.13]). Based on this, we can write Y(I h ), the variety 
generated by I h , as 

Y(I h ) = V(Ji) U V(J 2 ) U • • • U Y(J m ), 

where V( Ji) is the variety generated by the primary ideal Ji , and it is an irreducible component of the variety 
Y(I h ). We note that, given a generating set for the ideal I h , namely the Grobner basis { 3 },3 2 , • ■ • ,3s}, we 
can obtain a generating set for each primary ideal Ji in ([ 2 ]) by a computational algorithm, like the command 
“primaryDecomposition” in Macaulay2 package. In this way, we have the following procedure for finding 
the irreducible components of the PjRL for an irreducible transfer function G(s) = n(s)/d(s): 

1. Define p(s) = k d d(s) + k n n(s) and taking s = x+jy , obtain p{x+jy) = u(x, y , k d , k n ) +jv(x, y, k dl k n ); 

2. Obtain a Grobner basis for the ideal I = (u, v), and denote it by { 31 , 32 , ■ ■ ■, g s }', 

3. Let I h = ( 31 ,3 2 , ■ ■ ■ , 3 s), the homogenization of I, and obtain the primary decomposition of the ideal 
I h , as presented in Equation ([ 2 ]); 

4. The zeros of the generating set for each Ji in ([ 2 ]) is an irreducible variety, whose union for i = 1,2 ,,m 
makes up the PjRL. 


2.1 Examples 

In all examples below we used the Macalay2 software ([4]) to make the calculations and all polinomials 
are defined with coeficients in the field of rationals (that is in Q[xi, x 2 ,... ,x n \) so that we can get infinite 
precision in calculations. 

Example 2.1. Let be G(s ) = (s + l)/s 2 , whose RL plot is shown in Figure [ 2 J and 

p(s) = k d s 2 + k n (s + 1). 

Defining u = Re{p(ir + jy)} and v = Im{p(a; + jy)} we have: 

u(x, y, k d , k n ) = k d (x 2 - y 2 ) + k n (x + 1), v(x, y, k d , k n ) = 2 k d xy + k n y. 

Now we compute the Grobner basis for the ideal (u, v) using the graded reversed lexicographic order with 
x > y > k d > k n and obtain {31,32,33,34}, where: 


9 i(x,y,k d ,k n ) 

g 2 {x,y,k d ,k n ) 

g3(x,y,k d ,k n ) 

9 i{x, y, k d , k n ) 


2 xyk d + yk n (= r) 
x 2 k d - y 2 k d + xk n + k n (= q) 
x 2 yk n + y 3 k n + 2 xyk n 
2 y 3 k d - xyk n - 2yk n 


Homogenizing of the polynomials 3 ^, using the procedure indicated in the Introduction we obtain: 

31 = z 3 g 1 (x/z,y/z,k d /z,k n /z) = 2xyk d + yzk n 

32 = z 3 g 2 (x/z,y/z 1 k d /z, k n /z) = x 2 k d - y 2 k d + xzk n + z 2 k n 

33 = z A g 3 (x/z, y/z , k d /z , k n /z) = x 2 yk n + y 3 k n + 2xyzk n 

34 = z^g^x/z, y/z , k d /z, k n /z) = 2y 3 k d - xyzk n - 2yz 2 k n 
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Now we compute the primary decomposition for the ideal I h = (g 3 , g%, g 3 , g%) to obtain I h = J\ n J 2 D J 3 , 
where: 


■h 

J2 

J 3 


( x 2 +y 2 + 2 xz, 2 xkd + zk n ) 
(k r i- k n ) 


(3) 

(4) 

(5) 


Using the fact that kd and k n belong to the set of reals and that they can’t be both simultaneously zero, we 
have that 1 £ J 3 (suppose k n 7 ^ 0, so (1 /k n ) x k n = 1 € J 3 ) and then J 3 = R[x,y,z,kd,k n ] can be deleted 
from the primary decomposition of I h , that is, I h = J\ fl J 2 . Then we have that Y (I h ) = V(Ji) U V(J 2 ), 
where V( Ji) is defined by y = 0 and V( J 2 ) is defined by x 2 + y 2 + 2 xz = 0. To analyze these varieties in the 
affine XY plane we set z = 1 and we obtain the components of the plot shown in Figure [2j that is the line 
y = 0 and the circle (x + l ) 2 + y 2 = 1, as desired. Also, from the ideals J\ and J 2 above we can obtain the 
parametrization of these curves as well as the initial and terminal points of the PjRL, as was done by the 
author in [Tj: 

Variety V(Ji): Defined by the ideal Ji, as shown in Equation © 

• Initial Points: kd = 1 and k n = 0. From ©, we get y = 0 and x 2 = 0 or x = 0. Therefore the 
initial point for V( Ji) is (0 : 0 : 1) or (0, 0) in affine plane XY. 

• Terminal points: kd = 0 and k n = 1. We get from ©, y = 0 and xz + z 2 = 0. Then we have (a) 
z = 0 and x = 1, which is the point at infinity (1:0:0) (horizontal lines) and (b) z = 1 which 
implies x = — 1 and the point is (—1 : 0 : 1) or (—1,0) in affine plane XY. 

• Intermediary Points: kd = 1 and k n = A 7 ^ 0. Again from ©, we have y = 0 and x 2 +xz\ + z 2 X = 
0 , and we must have z = 1 (z = 0 would imply x = 0 what is impossible); so, all intermediary 
points are at finite position and is given by y = 0 and x 2 +xX+X = 0, A 7 ^ 0. Then x = ~ a±v/ 2 a2 ~ — ; 
which give us the RL over y = 0 line. 

Variety V(J 2 ): Defined by the ideal J 2 , as shown in Equation Q 

• Initial Points: kd = 1 and k n = 0. From ©, we get x 2 + y 2 + 2 xz = 0 and 2x = 0 or x = 0. Then 
we have x 2 + y 2 = 0 or y = 0. Therefore the initial point for V( J 2 ) is (0 : 0 : 1) or (0,0) in affine 
plane XY. 

• Terminal points: kd = 0 and k n = 1. We get from ©, x 2 + y 2 + 2 xz = 0 and z = 0. Then we 
get x 2 + y 2 = 0 what implies x = y = z = 0 which is not allowed, so this variety has no terminal 
points. 

• Intermediary Points: kd = 1 and k n = A / 0. Again from ©, we have x 2 + y 2 + 2 xz = 0 and 
2x + zX = 0, and we must have z = 1 (z = 0 would imply x = y = 0 what is impossible); so, all 
intermediary points are at finite position and is given by x 2 + y 2 + 2x = 0 and 2x + A = 0, A 7 ^ 0, 
which is the parametrized equation of the circle as shown in RL plot. 

Now, since the complete PjRL is V(Ji) U V(J 2 ), we have: 

• Initial points: (0:0:1) (from V(Ji)) plus (0:0:1) (from V(J 2 )); so we have a duplicate point at 
(0 : 0 : 1) or at (0, 0) in affine plane XY. 

• Terminal Points: {(1 : 0 : 0), (—1 : 0 : 1)}, only from V(Ji) 


Intermediary points: {(a; : 0 : l),x = ^ 0 } from V(Ji) plus {(x : y : l),x = —A/2, A 7 ^ 0 

and x 2 + y 2 + 2x = 0} from V(J 2 ). 


Example 2.2. We now consider a modification of Example 

RL plot is shown in Figure [3| Also in this case we see that 
“weird” curve shown in red. In this case: 


2.1 


above by defining G(s) = 


s + 1 


-, whose 


s 2 (s + 4)' 

the RL is the union of the line y = 0 and the 


and we have 


p(s) = k d (s 3 + 4s 2 ) + k n (s + 1) 


u = Re{p(x + jy)} = kd{x 3 — 3xy 2 + 4x 2 — 4y 2 ) + k n (x + 1) 
v = Tm{p(x + jy)} = k d {-y 3 + 3 x 2 y + 8 xy) + k n y. 
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The Grobner basis for the ideal (u, v) is {gi,g 2 , < 73 , < 74 , 35 }, and the generating set for I h is {g ^, <7 2 , <73 , < 74 , <75 }, 
where: 


9i = 3x 2 yk d - y 3 k d + 8a :yzk d + yz 2 k n 
g% = x 3 k d — 3xy 2 k d + Ax 2 zk d — Ay 2 zk d + xz 2 k„ + z 3 k n 
g% = 2 x 3 yk n + 2 xy 3 k n + 7x 2 yzk n + 3 y 3 zk n + 8xyz 2 k n 
g'l = 2Axy 3 k d + 32y 3 zk d + 32xyz 2 k d - 6xyz 2 k n - hyz 3 k n 

g§ = 24 y 5 k d + 160 y 3 z 2 k d — 128xyz 3 k d — 18a : 2 yz 2 k n — 24 y 3 z 2 k n — 39 xyz 3 k n — 52 yz A k n 
and computing the primary decomposition for I h we obtain I h = Ji D J2 D J3, where 


Ji = (y, x 3 k d + Ax 2 zk d + xz 2 k n + z 3 k n ) 

J 2 = {2x 3 + 2xy 2 + 7x 2 z + 2>y 2 z + 8xz 2 , 3x 2 k d — y 2 k d + 8xzk d + z 2 k n ) 

J3 — {k d i k n ^j 


( 6 ) 

(7) 

( 8 ) 


Once more, J3 can be removed from the intersection, so Y (I h ) = V(Ji) UV(J 2 ), wh ere V(Ji) is defined by 
y = 0 and V(J 2 ) is defined by 2a ; 3 + 2a;j / 2 + 7 a : 2 2 + 3?/ 2 2 + 8 a: 2: 2 = 0. As in Example 2d]above, to analyze these 


varieties in the affine XY plane we set z = 1 and we obtain the components of the plot shown in Figure [3j 
that is the line y = 0 and the curve plotted in red whose equation is 2a; 3 + 2xy 2 + 7x 2 + 3 y 2 + 8x = 0. We 
can also plot these varieties over the projective plane (using Gnomonic projection) as was shown in jlj, as 
well as plot them over other components of the projective plane also, as the affine plane ZY. 


3 Algebraic Dual Root-Locus — ADRL 

Duality is a fundamental concept in projective algebraic geometry. In fact, it is a basic property of the real 
projective plane MP 2 that a “point” with (nonzero) homogeneous coordinate (a : b : c) can be associated 
to a “line” with equation ax + by + cz = 0 in and vice-versa. This kind of duality can be extended from 
a projective line to a projective plane curve (W) defined by f(x,y,z ) = 0, where f(x,y,z) £ M[a;,y,z] is a 
homogeneous polynomial. By the natural duality between lines and points in RP 2 , each tangent line to the 
curve can be associated to a point with homogeneous coordinate, for instance, (u : v : w), and the main result 
is that this set of points is also the solution to some equation g(u,v,w) = 0, where g(u,v,w ) £ M[it, v,w\ is 
a homogeneous polynomial, that represents a curve (W*) over RP 2 (in fact over the dual of RP 2 , which it 
is itself). Therefore W* is denominated the dual curve of W. Interestingly, we also have that if we take the 
dual of the dual of a curve we restore back the original curve, that is (W*)* = W (see [8)). Mathematically, 
the dual curve of f(x,y,z) = 0 is the set of points (u : v : w) = (df/dx : df/dy : df/dz) in PR 2 , or 
u = A (df/dx), v = \(df/dy) and w = A (df/dz), for some A / 0. To find the curve which these points 
belongs to, we can restate the problem as the one of eliminating x , y, z and A from the set of equations 
f(x, y, z) = 0, u — A (df/dx) = 0, v — \(df/dy) = 0 and w — A(< df/dz) = 0, which, in turn, can be solved by 
finding a Grobner basis for the ideal 

T / f ( \ \ d f \ d f \ d f\ (Q\ 

I = (*, y,z),u- X-v-X-,w-X-J ( 9 ) 
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In our context, we are primarily interested in obtaining the dual curve for each irreducible component of 
the RL, as computed in Section [2] above, and collate them to construct a new RL plot that we denominate 
“Algebraic Dual Root-Locus” (ADRL). So, we will calculate de ideal defined in ([9]) for the examples presented 
in Section [3 


3.1 Examples 

Example 3.1. Let be G(s) = (s + l)/s 2 , whose RL plot is shown in Figure[2j and we found in Example 
that the PjRL can be represented by the set of ideals in Equations (|3|4[), that is 


2.1 


Ji = (y, x 2 kd + xzk n + z 2 k n ) 

Ji = (x 2 + y 2 + 2 xz, 2xkd + zk n ) 


( 10 ) 

( 11 ) 


To find the dual curve of V( Ji), we have to calculate the dual curve of /(x, y, z) = x 2 kd + xzk n + z 2 k n , 
and therefore the ideal in Equation (|9| for this curve is: 

I = ( x 2 kd + xzk n + z 2 k n , u — X(2xkd + zk n ),v, w — X(xk n + 2 zk n )) 

Finding a Grobner basis for this ideal with “Lex” monomial ordering with x>y>z>kd>k n >X>u> 
v > w we eliminate x, y, z, X and obtain the curve: 

g(u, v , w) = kdW 2 + k n u 2 — k n uw 

So we can define the “dual ideal” for J\ as: 

jf = (v, kdW 2 + k n u 2 — k n uw) 


( 12 ) 

2 1 2 xz, and 


and the dual curve for Y(Ji) is Y(Jf). 

To find the dual curve of Y( J 2 ), we have to calculate the dual curve of f(x, y, z) = x 2 + y 
therefore the ideal in Equation ([ 9 ]) for this curve is: 

I = (x 2 + y 2 + 2 xz, u — A(2x + 2 z),v — A(2 y),w — A(2x)) 

Finding a Grobner basis for this ideal, to eliminate x, y, z, A, we obtain 

h(u, v, w) = v 2 + 2uw — w 2 

To find the parametrization for h above we use the ideal as shown below 

I = (2 xkd + zk ni u — X(2x + 2z), v — A(2y), w — A(2x)) 

and finding a Grobner basis for this ideal, to eliminate x, y, z, A, we obtain the (parametrized) curve 
hi(u, x, w) = k n u + (2 kd — k n )w and then we can define the again “dual ideal” for J 2 as: 


J -2 = ( y2 + 2 ilw — w 2 , k n u + (2 kd — k n )w) 


(13) 


Now we can define the Algebraic Dual Root-Locus for G(s ) = (s + l)/s 2 as the union of the varieties V(Jf) 
and V(J|) where Jf and are defined in (12) and (13), respectively. Below we analyse each variety in 
order to obtain the complete plot for the ADRL: 


Variety V(J^): Defined by the ideal Jf, as shown in Equation (12) 


• Initial Points: kd = 1 and k n = 0. From (12), we get v = 0 and w 2 = 0 or w = 0. Therefore the 
initial point for Y(Ji) is (1 : 0 : 0), a point at infinity, or the intersection of horizontal lines in 
affine UV plane. 

• Terminal points: kd = 0 and k n = 1. We get from (12), v = 0 and u 2 — uw = 0. Then we have 
(a) u = 0 and w = 1, which is the point (0:0: 1), or the point (0,0) in affine plane UV and (b) 
u = w = 1 which is the point (1 : 0 : 1) or (1, 0) in affine UV plane. 

• Intermediary Points: kd = 1 and k n = A / 0. Again from (12), we have v = 0 and Xu 2 — 
Xuw + w 2 = 0, and we must have w = 1 (w = 0 would imply u = 0 what is impossible); so, all 
intermediary points are at finite position (w = 1) and is given by v = 0 and u 2 — u + 1/A = 0, 

A ^ 0. Then u = which give us the ADRL over v = 0 line in affine UV plane. 


Variety Defined by the ideal J.J, as shown in Equation (13) 
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Figure 4: Algebraic Dual Root-Locus for G(s) = (s + l)/s 2 


Initial Points: k d = 1 and k n = 0. From (131, we get v 2 + 2uw — w 2 = 0 and 2w = 0 or w = 0. 
Then we have v 2 = 0 or v = 0. Therefore the initial point for ¥(^ 2 ) is (1 : 0 : 0), a point at 
infinity (intersection of horizontal lines) in UV plane. 


Terminal points: k d = 0 and k n = 1. We get from (13), v 2 + 2 uw — w 2 = 0 and u — w = 0 or 
u = w. Then we get v 2 + u 2 = 0 what implies u = v = w = 0 which is not allowed, so this variety 
has no terminal points. 


Intermediary Points: k d = 1 and k n = A / 0. Again from (13), we have v 2 + 2 uw — w 2 = 0 and 
Xu + (2 — A )w = 0, and we must have w = 1 (w = 0 would imply u = v = 0 what is impossible); 
so, all intermediary points are at finite position and are given by v 2 + 2u— 1 = 0 and it = 1 — 2/A, 
A 7 ^ 0, which is the parametrized equation of the parabola as shown in ADRL plot. 


Now, since the complete ADRL is V(jf) UV(J 2 d ), we have: 

• Initial points: (1:0:0) (from V(J/)) plus (1:0:0) (from V/J^)); so we have a duplicate point at 
(1 : 0 : 0) or at infinity in affine plane UV. 


• Terminal Points: {(0 : 0 : 1), (1 : 0 : 1)} only from V(j/); or (0.0) and (1,0) in UV plane 


• Intermediary points: {(it : 0 : l),u = (1 ± y/l — 4/A)/2, A 7 ^ 0} from V(J-f) plus {(u : v : l),u = 
1 — 2/A, A 7 ^ 0 and v 2 + 2u — 1 = 0} from V(T|), in UV plane. 


The ADRL plot for G(s) = (s + l)/s 2 is shown in Figure |4j 

It is important to note that if we take the duals of Jf and J 2 , defined in Equations (12) and (13), 
respectively, we get back, the ideals J\ and J 2 , as defined in Equations (10) and respectively. The 
procedure for doing that is eliminating u,v, w and A in the ideals Ii, I 2 and I 3 (defined below), using Lex 
monomial ordering with u>v>w>x>y>z>k d >k n >\. 


h = ( k d w 2 + k n u 2 — k n uw, x — X(2uk n — k n w), y,z — A(2 k d w — k n u)) 

1 2 = (u 2 A 2uw — ic 2 , x — A(2ic), y — A(2u), z — \{2u — 2ic)) 

1 3 = (k n u + (2k d — k n )w,x — X(2w)y — X(2v),z — X(2u — 2w)) 


Example 3.2. We now consider G(s) = 


s + 1 


represented by the ideals in Equations (6|7): 


(s + 4) 


whose RL plot is shown in Figure jsj and whose PjRL is 


Ji = (y, x 3 k d + Ax 2 zk d + xz 2 k n + z 3 k n ) 

J 2 = (2x 3 + 2xy 2 + 7x 2 z + 3y 2 z + 8xz 2 ,3x 2 k d —y 2 k d + 8xzk d +z 2 k n ) 

Repeating the reasoning used in Example |3.1| above, we obtain the “duals” ideals: 

J-f = (v, 4:k d uw 2 — k d w 3 + k n u 3 — k n u 2 w) 

J 2 = (f(u,v,w),g{u,v,w)) 
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(14) 

(15) 














where / and g are the “astonishing” polynomial^] 

f(u, v , w) = 216m 5 w + 144m 4 m 1 2 — 621m 4 u> 2 — 912u 3 v 2 w + 720m 3 w 3 — 352m 2 m 4 + 7\8u 2 v 2 w 2 — 

424m 2 m> 4 — 232 uv A w — 176mm 2 m; 3 + 128mm; 5 — 240m 6 + 107 v 4 w 2 — 8 v 2 w 4 — 16m; 6 
g{u, v, w) = 294912/cJm 3 m; - 327680/cJmV - 798720kju 2 w 2 + 1671168 k\uv 2 w + bUOOOk^uw 3 - 
1048576fc^M 4 - 675840fc> 2 M; 2 - 96000fc> 4 + 73728k 3 ,k n v 4 - 448512fc^/c n M 3 M; + 

430080A; d k n u 2 v 2 + 706304k d k n u 2 w 2 — 788480k 3 d k n uv 2 w — 403200 k d k n uw 3 + 
131072^fc„M 4 + 300288k d k n v 2 w 2 + 75600 k 3 d k n w 4 - 46656/^fc 2 M 4 + 174816A;^ 2 m 3 m; - 
68928k d k 2 u 2 v 2 - 224292k 2 d k 2 n u 2 w 2 + I09920k 2 d k 2 n uv 2 w + l\9040k 2 d k 2 n uw 3 - 
6144fc^ 2 M 4 - 41364/c^fc 2 M 2 tM 2 - 22320 k 2 d k 2 n w 4 + 8208k d klu 4 - 26172/c d fc 3 M 3 u; + 

3088fc dfc 3 M 2 v 2 + 30636 k d k^u 2 w 2 — 4764k d k^uv 2 w — 15616k d k^uw 3 + 

128 k d k 3 v 4 + I788k d k^v 2 w 2 + 2928 k d k 3 w 4 — 441k 4 u 4 + 1344fc 4 M 3 u> — 

42k 4 u 2 v 2 — 1528fc 4 M 2 M) 2 + 64k 4 uv 2 w + 768k 4 uw 3 — k 4 v 4 — 24fc 4 u 2 u; 2 — 144 k 4 w 4 

Remembering that the dual curve for 2a; 3 + 2 xy 2 + 7x 2 z + 3 y 2 z + 8xz 2 = 0 is f(u,v,w ) = 0 and that 
g(u,v,w) = 0 is used just for obtaining a parametrization. The plot for the correspondig ADRL, obtained 
by the Scilab package ([5]) is shown in Figure [5j 


4 Conclusions 

We have presented in this paper a procedure for isolating the planes curves that makes up the root-locus 
plot for an irreducible transfer function. This procedure can be easily implemented in a computational 
algebra software package. We also showed how to compute the dual curve (in projective algebraic geometry 
sense) to each plane curve and join them to compose what we denominated “Algebraic Dual Root Locus” 
(ADRL). Some examples were worked out in order to show the effectiveness of the procedure. We intend to 
investigatethe properties of the ADRL more deeply in future works. 

1 In fact, a well known result in projective algebraic geometry states that if a curve has degree d (and no singularities) its 

dual has degree d{d — 1) [8j pp. 173]; in our case the polynomial in the ideal J 2 has degree 3, so its dual has degree 6. 
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